Branching-annihilating random walks in one dimension: Some exact results 



K. Mussawisade 1 , J. E. Santos 2 and G. M. Schutz 1 

1 Institut fur Festkorperforschung, Forschungszentrurn Jiilich, 52425 Jiilich, Germany 
2 Department of Theoretical Physics, University of Oxford, 1 Keble Road, Oxford, 0X1 3NP, United Kingdom 

(February 1, 2008) 

We derive a self-duality relation for a one-dimensional model of branching and annihilating random 
walkers with an even number of offsprings. With the duality relation and by deriving exact results 
in some limiting cases involving fast diffusion we obtain new information on the location and nature 
of the phase transition line between an active stationary state (non-zero density) and an absorbing 
state (extinction of all particles), thus clarifying some so far open problems. In these limits the 
transition is mean-field-like, but on the active side of the phase transition line the fluctuation in the 
number of particles deviates from its mean-field value. We also show that well within the active 
region of the phase diagram a finite system approaches the absorbing state very slowly on a time 
scale which diverges exponentially in system size. 
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I. INTRODUCTION 

In a branching-annihilating random walk (BARW) 
particles hop on a lattice, annihilate pairwise on en- 
counter, but may also spontaneously create offsprings on 
the same or on nearest neighbour lattice sites. Such mod- 
els appear in a large variety of contexts, in particular 
in reaction-diffusion mechanisms and in non-equilibrium 
spin relaxation. A generic feature of these processes is a 
transition as a function of the annihilation and branch- 
ing rates between an active stationary state with non-zero 
particle density and an absorbing, inactive state in which 
all particles are extinct. Numerical results gained from 
a large variety of systems suggest that the transition in 
models with a single (or an odd number) of offsprings fall 
into the universality class of directed percolation (DP) 
0, whereas models with an even number of offsprings 
belong to a distinct parity-conserving (PC) universality 
class d|| . A coherent picture of this scenario is provided 
from an renormalization point of view Q . 

In this paper we use exact methods to derive a self- 
duality relation and to address some open questions for 
limiting cases of the BARW model of Ref. Q which is 
a model for spin relaxation dynamics far from thermal 
equilibrium. One considers Ising spins in one dimen- 
sion with generalized zero-temperature Glauber dynam- 
ics [^J, but with an independent coupling to an infinite- 
temperature heat bath which allows for Kawasaki spin- 
exchange events Q with rate a/2. This spin- flip process 
can be visualized in the following way: 



t i T -» T T T and | T I -» I I 1 with rateA 
T T I - T I I and 1 1 T - I T T with rate D/2 (1) 
TI^IT and|T^T4 with rate a/2 

By identifying a domain wall (f j, or J, |) with a particle 
of type A on the dual lattice and two parallel spins with 
a vacancy Jt| this process becomes a BARW with rates 

A A —> with rate A 

% A ^ A% with rate D/2 (2) 
A ^ AAA with rate a/2 
% A A ^ A A with rate a/2 

Without branching (i.e. spin exchange) the system 
evolves into the single absorbing state with no particles 
at all. In spin language this is the totally ferromagnetic 
state with all spins up or all spins down. In the pres- 
ence of the branching process an intricate competition 
between the zero-temperature ordering process (parti- 
cle annihilation) and the disordering high-temperature 
branching process sets in. The result is a non-trivial 
phase diagram as a function of the system parameters. 
Starting from, say, a random initial state with an even 
number of particles the system evolves ultimately into the 
inactive empty lattice for dominant ordering dynamics, 
whereas it remains in an active state with finite density if 
the disordering branching process dominates. Numerical 
evidence suggests the phase transition to belong to the 
PC universality class ||. 

We stress that these results are supposed to be valid 
only in the thermodynamic limit. In any finite system 
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the unique stationary state is the absorbing inactive state 
(unless A = 0) because for A > there is always a small 
probability of reaching this state from which the system 
cannot escape any more. However, intuitively, one ex- 
pects the approach to this state to occur on a time scale 
Tact ~ Q L which is exponentially large in system size L 
if parameters are chosen to represent the active phase of 
the thermodynamic limit. In the absorbing phase, both 
exact analytical results for A = D || and renormalization 
group results on diffusion-limited annihilation {a — 0) |IJ 
show that the approach to extinction is algebraic for the 
infinite system. For a finite system one can infer from 
these results a crossover time scale T a bs ~ L 2 to expo- 
nential decay of the particle density. 

Here we aim at obtaining information on the form of 
the phase transition line and on the dynamical and sta- 
tionary behaviour of the system in various limiting cases 
involving fast diffusion of particles (or spins in terms of 
the spin-relaxation model). In Sec. [n] we derive a du- 
ality relation which maps the phase diagram onto itself 
in a non-trivial way. This duality is different from the 
domain-wall duality which maps the spin-relaxation dy- 
namics to the BARW. On a self-dual line running across 
the phase diagram we obtain a relation between the 
(time-dependent) density expectation value for a half- 
filled random initial state and the survival probability of 
two single particles in an initially otherwise empty sys- 
tem. 



In Sections [II and IV we adopt the strategy of con- 
sidering separately the fluctuations in the total number 
of particles from the spatial correlations within a con- 
figuration of a given fixed number of particles. By sep- 
arating the hopping time scale from the branching and 
annihilation time scales one can then gain insight in the 
behaviour of the system in the absence of spatial cor- 
relations. In Sec. Ill we study the system i n the fast- 
diffusion limit D -► oo of the BARW (§) (|IIIA|) and 
in its spin-relaxation formulation (^) in the dual limit 
of infinite spin-exchange rate a. In these limiting cases 
all spatial correlations are washed out and one expects 
the PC transition to change into some other, mean-field- 
type phase transition. However, in contrast to a tradi- 
tional mean-field approach, our treatment keeps track of 
the exact fluctuations of the total number of particles (or 
spins respectively). Our treatment is not an approxima- 
tion, but yields rigorous results in these limiting cases 
for which we calcu late the stationary density and density 
fluctu ations ( III A ) and fluctuations of the magnetization 
( III B ). Thus we are able to analyze to which extent the 
system deviates from mean-field behaviour and we also 
identify the exact phase transition point. In Sec. IV we 



investigate by similar means the dynamical behaviour of 
the finite system in the active region of the phase dia- 
gram. We show that for fast diffusion the relaxation to 
the absorbing state in a finite system is indeed exponen- 
tially slow, thus confirming the intuitive argument for the 
signature of the active region in a finite system. 
In Section [v| we conclude with some final remarks. 



II. DUALITY RELATIONS 

We define the BARW in terms of a master equation 
for the probability P(rj; t) of finding, at time i, a con- 
figuration rj of particles on a lattice of L sites. Here 
V = {vO-)iV(2)> • • • >il(L)} where r](x) = 0,1 are the 
integer-valued particle occupation numbers at site x. For 
definiteness we assume L to be even. Using standard 
techniques [ ft0| we express the time evolution given by 
the master equation in terms of a Hamiltonian H of 
a one-dimensional quantum spin system (for a review 
see 0). The idea is to represent each of the possible 
particle configurations 77 by a column vector 1 77 ) which 
together with the transposed vectors ( 77 1 form an or- 
thonormal basis of a vector space X = (C 2 )® L . One 
represents the probability distribution by a state vector 

I P{t) ) = J2riex *)l V) an d writes the master equa- 
tion in the form 



-P( m t) = -(r,\H\P(ty : 



(3) 



where the off-diagonal matrix elements of H are the (neg- 
ative) transition rates between states and the diagonal 
entries are the inverse of the exponentially distributed 
life times of the states. A state at time t' = to + t is 
given in terms of an initial state at time to by 

\P(t + t))= C - Ht \P(t )). (4) 

The expectation value Pk(t) = (s\rik \ P{t) ) for the den- 
sity at site k is given by the projection operator rik which 
has value 1 if there is a particle at site k and otherwise. 
The summation vector ( s \ = X^ex ( V I performs the av- 
erage over all possible final states of the stochastic time 
evolution. Below an initial distribution with N parti- 
cles placed on sites k%,. ..,kjy is denoted by the column 
vectors | k\, . . . , kjsr )■ The empty lattice is represented 
by the vector | ) . The uncorrelated product distribu- 
tion where on each lattice site the probability of finding 
a particle is equal to 1/2, is given in terms of the trans- 
posed of the summation vector as 1 1/2) = ( s | /(2 L ). 

To obtain the Hamiltonian for the time evolution of 
the BARW (|^) we note that we can represent any two- 
state particle system as a spin system by identifying a 
particle (vacancy) on site k with a spin-up (down) state 
on this site. This allows for a representation of H in 
terms of Pauli matrices where n& = (1 — crf)/2 projects 
on states with a particle on site k and Vk = 1 — rik is the 
projector on vacancies. The off-diagonal matrices s k = 
(a^ ± iaT)/2 create (s7) and annihilate (s k ) particles. 
We stress that in the present context the "spins" are 
just convenient labels for particle occupancies which are 
conceptually entirely unrelated to the spins of the spin 
relaxation model ([!]) which is treated below. Using this 
pseudospin formalism one finds 



H 



rikVk+i + VkUk+i - s k s k+1 



Si, S 



+ 



k °k+l> 



+2\{n k n k+ i - s+ s k+i) + a (! ~ °k<4+i) n k} ■ (5) 
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Each part of this stochastic Hamiltonian represents one 
of the elementary processes (§) and we may write 

H(D, A, a) = DH SEP + \H RSA + aH BARW . (6) 

Here H SEP represents hopping of hard-core particles, 
i.e. the symmetric exclusion process |l2] ], the pair- 
annihilation process encoded in H RSA corresponds to 
random-sequential adsorption |Tj| and H BARW describes 
a special equilibrium branching-annihilating random 
walk where pair annihilation requires the presence of 
another particle and hopping occurs only in pairs (see 
(^). For A = D and a = this process reduces to 
the exactly solvable process of diffusion-limited pair an- 
nihilation (DLPA) |l4|| . The time evolution conserves 
particle number modulo 2. Here we work only on the 
even subspace defined by the projector (1 + Q)/2 where 
Q = (— 1)^ = rifc^fc- The projection on the even sec- 
tor of the uncorrelated initial state with a density 1 /2 is 
given by the vector | l/2) even = (1/2) L - 1 | s ) even . 

Within the same framework the stochastic Hamilto- 
nian for the spin-flip process ([!]) can be written in terms 
of Pauli matrices as follows 



k 

+a(l-<7f<7f +1 )(l-C7^ +1 )] 



(7) 



with the generalized Glauber spin flip rates encoded 
in w k (D,X) = (2 - al_ x al - <r z k a* k+1 )(D + A + (A - 
£ > )c r fe_i cr fe + i)/2. For this process the spins represent 
the actual spin configurations of the spin-relaxation pro- 
cess. We note that the usual zero-temperature Glauber 
dynamics - equivalent to DLPA in particle language - 
correspond to A = D, a = 0. For this model the do- 
main wall correspondence between the BARW and the 
spin-flip process can be rigorously derived as a similarity 
transformation on the level of the quantum Hamiltonian 
description. There exists a transformation B such that 
H SF = BHB- 1 fT|. The generalization A = D, a > 
corresponds to the exactly solvable process introduced in 
Ref. §. 

Consider now the transformation T>± which is, for the 
even particle sector, defined by 



£>+ = 7l72 • ■ ■ 72L-1 



where 



72fc „! = - [(1 + i)a% - (1 - »)] 

72fc = - [(1 + *K<7fe+l - (1 - i)] 



(8) 

(9) 
(10) 



and defined by T>_ = — r D + a' J [ for the odd particle sector 
|^6|,0. V± is unitary and transforms Pauli matrices as 
follows: 



' D ± o-%a% +1 V ± 



u k 



k = L 



^k+i k + L 
±Qo x L o\ k = L 



(11) 
(12) 



In Ref. |lq| it was observed that this transformation 
maps H DI ^ r ^, obtained from H by setting A — D and 
a = 0, onto its transposed, H DLPA = V{H DLPA ) T V~ 1 
and thus generates a set of relations between various ex- 
pectation values [^9|. Here we go further and apply this 
transformation to the Hamiltonian H = H(D, A, a) 
and transpose the operator which results from the trans- 
formation. Using 0), (H) we find 



H 



H SEP_ 


^ j^BARW 


BARW _ 


+ H SEP , 


tjRSA _ 


-> H RSA + 



H 



SEP 



H 



BARW 



and hence the relations 

H = XH RSA + (A + a)H SEP + (D - \)H BARW . (13) 

which has the same form as the original Hamiltonian (|TJ) , 
but with rates 



A = A, 
D = A + a, 
a = D - A. 



(14) 



The transformation ( |13[ ) is a duality transformation, we 
obtain the identity transformation if we apply the trans- 
formation twice. 

On the mean-field phase transition line A = D [^o| the 
system is exactly solvable [§j and belongs in its entirety 
to the inactive phase. Hence the whole region A > D is 
in the inactive phase. Duality maps the interesting re- 
gion A < D which contains the phase transition line non- 
trivially onto itself. Thus duality can be used to relate 
physical quantities at different points of the parameter 
space. This region contains a self-dual line 



D = A 



(15) 



in which every point maps onto itself. In the notation of 
Ref. @ D = Prw = T-^l-S), A = Pan = r- 1 (l + <5) and 
a = 2p ex = 2(1 — 2r _1 ) is normalized such that p ex + 
Prw +Pan = 1< With a parametrization in terms of 8 and 
p ex the dual rates are given by S — —2p ex /[l+p ex + 6(l — 
p ex )} andp~ ex = -S(l-p ex )/[S(l-p ex )+2(l+p ex )]. The 
self-dual line is given by the relation i5 = —2p ex / (l—p ex ). 

The duality transformation not only maps the 
phase-diagram onto itself, but also generates rela- 
tions between time-dependent expectation values. Con- 
sider the expectation value of the density p k (t) — 
(s\ even rike~ Ht \l/2) even , where \ l/2) even is a random 
initial state with density 1/2, projected over the even 
sector. This expectation value is defined at the point 
(D, A, a) of the parameter space of the Hamiltonian. It 
is straightforward to verify the relations 



{l/2\ em 



' en = -i(i-l) L - l \ 
l V = i(-i-l) L - 1 . 



o; 

0|/2 



L-l 



(16) 
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So if we use these rules of transformation and the rules 



of transformation for the Pauli matrices, given by (11), 
(|l2|), we can write the expectation value for the density 
in the even sector the form 



S Ke- at |l/2r en = ^(0|e- at (l 



ok 



'fe-i°fe/ 



-Ht\ 



fc, fe + 1 



(17) 



where we have used (13) and the fact that the expecta- 
tion value for the density p%[t) is a real number. The 
transformed initial state is a superposition of the steady 
state (the empty lattice) and the two-particle state with 
particles at sites k, k+ 1. The quantity in the right hand 
side of equation ( p7| ) is one-half times the probability 
that the state with two particles initially placed at sites 
k and k + 1 has not decayed at time t to the empty state, 
measured with the transformed rates A, D, a (|14|). 

This is a specific result for the time-dependent den- 
sity starting from a random initial state with density 
1/2. More general transformation properties of time- 
dependent correlation functions can be obtained follow- 
ing the strategy of Ref. We conclude this section 
by pointing out that analogous enantiodromy relations 
can be derived for a discrete-time version of the pro- 
cess which corresponds to a sublattice parallel updating 
scheme rather than the random sequential updating rep- 
resented by the stochastic Hamiltonians (^|), (0). Such 
a parallel updating scheme (which we expect to retain 
all the universal features of the model) consists of four 
steps. In a first updating sweep, update all spins on the 
even sublattice in parallel according to the generalized 
Glauber rules, but with the rates A, D now taken as ac- 
tual probabilities. In a second step one updates the odd 
sublattice. In a third step one applies a sublattice parallel 
pair-updating scheme to implement the Kawasaki spin- 
exchange process with probability a/2: In a first round 
one divides the lattice into even/odd pairs (2k, 2k + 1) 
and exchanges spins within each pair with probability 
a/2. Finally one updates with the odd/even pairs. This 
completes a full updating cycle. The stochastic time evo- 
lution of this process may be written in terms of the 
transfer matrix 

T SF = T D K dd (a) T L K ven (a) T Q G d(i (A, D) Tg en {\, D) (18) 

where 

L/2 

Tg en (X,D) = II [1 - (1 - o%)w k {D,\)/4] (19) 

k=l 

and an analogous expression for T^ dd . The spin exchange 
transfer matrix T^ dd T^ en is the well-known transfer ma- 
trix of the six-vertex model J2l[] defined on a diagonal 
square lattice. One has 

L/2 

T* m = ft [1 - a(l - ^Vffc+iX 1 - °f**f*+i)/4] • 



k=l 



The transfer matrix for the related BARW model can 
be obtained by applying the similarity transformation 
B of Ref. [^5|. One can then derive a duality relation 
in the way described above. The transformed process 
has the same elementary transitions, but with a different 
updating sequence. 



III. PHASE TRANSITION FOR FAST 
DIFFUSION 

It is intuitively clear that the PC phase transition in 
the system originates in the complicated structure of the 
density correlations which are build up by the compet- 
ing processes of branching and annihilation. For a better 
understanding, consider first A = 0. This reduced pro- 
cess includes (besides diffusion) branching A — ► 3A and 
conditional pair annihilation 3A — > A which both require 
the presence of a surviving particles to take place. As a 
result, there are two stationary distributions: the empty 
lattice, and the random distribution where each particle 
configuration is equally likely. Since there is no transition 
channel from the occupied lattice to the empty lattice, 
the system is in the active phase. 

On the other hand, for a — the system is in the 
absorbing phase, the only stationary distribution is the 
empty lattice. We conclude that the unconstrained pair 
annihilation process 2A — > with rate A is responsible 
for the phase transition to take place. 

This scenario is captured in a simple mean-field ap- 
proach. The exact equations of motion for the expected 
particle number (N(t)) reads 



1( N(t) >=£>< > - 2(a + A)( ni (t)n i+1 (t) )] . 

i 

(21) 



Replacing the correlators by the product of the density 
p(t) = ( n, (t) ) yields the mean- field equation 

j t {N(t)) = a(N(t))- ^±^{N(t)f (22) 

with the stationary mean-field solution for the active 
phase 



2(q + A) 



(23) 



Since each lattice site can take only one particle and 
therefore rii = 0,1, one can use nf = m to show that 

the mean-field fluctuations A; nf = ( )* mf - ({ N )* n/ ) 2 
around the mean are given in terms of the density p* m ^ = 
(N)* mf /L by 



(24) 



(20) 



We conclude that the mean-field phase transition point 
is given by a/ A — 0, consistent with the considerations 
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above. By duality (jlj) we also recover the mean-field 
phase transition line of Ref . . 

There are several questions that we want to address 
in this context. The first is the form of the exact phase 
transition line if both a and A are very small compared to 
the diffusion rate D. The second question is the nature 
of the phase transition in this limit. If D 3> a, A the spa- 
tial correlations build up by the annihilation/branching 
process are wiped out very quickly by diffusive mixing, 
leaving a transition which we cannot expect to be a PC 
transition anymore. Finally, in the next section we study 
the crossover time scales on which a large, but finite sys- 
tem reaches the absorbing state. 



A. Phase transition in the BARW 

To tackle these questions we observe that for fast diffu- 
sion the process simplifies dramatically: In the absence of 
spatial correlations the state of the system is fully char- 
acterized by the total particle number N. For fixed N, 
each particle configuration occurs with equal probability 
Nl(L—N)\/(L)\ which is just the inverse of the number of 
possibilities of placing TV particles on a lattice of L sites. 
As a result, the dynamics reduce to a random walk on the 
integer set 0, 2, 4, ... , 2A", . . . , L of total particle number 
N = 2K. Thus we may represent the dynamics as a ran- 
dom walk on a one-dimensional lattice of L/2 + 1 sites, 
where the position of the random walker marks the num- 
ber of particles of the BARW process and 0, representing 
the empty lattice, is an absorbing point. It remains only 
to calculate the hopping rates rjv and In from site N to 
the right (N + 2) and left (N-2) respectively. The state 
of the system is then given by the solution of the master 
equation 



dt 



P N (t) = r N - 2 P N _ 2 (t) +e N+2 p N+2 (t) 

~(r N +£ N )P N (t) 



(25) 



for the probability of finding N particles in the sys- 
tem. The average particle number is given by ( N(t) ) = 

V v A7\.:/i f§. 

By counting the number of possibilities of finding two 
vacancies on neighbouring sites of an occupied site in a 
random state of N particles one readily finds 



a N(L — N)(L — N — 1) 

r N = 77" 



(L-l)(L-2) 



(26) 



as contribution from the branching process with rate a/2. 
An analogous consideration gives 



In = 



1 



aN{N - l)(N- 2) 2XN{N-1) 
(£-!)(£ -2) + (L-l) 



(27) 



as the contribution from the annihilation processes 3 A 
A and 2A — ► respectively. 



These rates represent a biased random walk which in 
the thermodynamic limit L — ► oo and for N fixed reduces 
to a directed random walk in positive direction with in- 
creasing hopping rate r^r = aN/2. Since for D — > oo 
the branching process is not diffusion-limited, the par- 
ticle number increases exponentially in time, ( N(t) ) = 
(N(0))e at . Thus for any a > the system is in the ac- 
tive phase, i.e., the phase transition is at a = which is 
consistent with the mean- field result (p3|). For a large, 
but finite system with a small initial number of particles 
one expects a slowing down of the exponential growth 
when a finite density is reached, i.e. on a time scale of 
the order In [L)/a. Ultimately, though, the finite system 
will reach, by a rare fluctuation in the number of parti- 
cles, the absorbing empty lattice. This second crossover 
time to absorption is discussed in the next section. 

To study the stationary behaviour d/(dt)P]sr(t) = of 
the system we rescale the lattice to unit length and ex- 
pand the r.h.s. of the master equation j25|) in a Taylor 
series in the lattice spacing 1/L. Setting x = N/L and 
keeping the leading order term yields the equation c = 
2(lx — r x )P* . The integrability condition f Q dxP* = Ion 
the stationary probability distribution requires the inte- 
gration constant c to vanish. The resulting equation has 
the solution P* = S(x), corresponding to the absorbing 
state. The only other integrable solution on the interval 
[0, 1] is the delta-function P* — 8(x - p*) with 



(28) 



A)' 



This gives the exact stationary density p* of the active 
phase which, not very surprisingly, coincides with the 
mean- field value (p3|). 

To determine whether the system in the infinite- 
diffusion limit actually is a mean-field system we inves- 
tigate the fluctuations around the mean J2jj|) . The mean 
field result ( p4| ) requires studying the fluctuations on a 
length scale of order y = y/X(x— p*). Keeping in the Tay- 
lor expansion of the master equation around x = p* all 
terms to this order gives the ordinary differential equa- 
tion 



dy v 



y 



2p*(l-p*) 



*\2 V 



(29) 



The solution of this equation is a Gaussian which gives 
the exact fluctuations in the particle number 



A* = 2/(1 -p*) 2 L 



(30) 



Except for A = (p* = 1/2) this expression is in dis- 
agreement with the mean- field result (f24j), indicating a 
non-trivial effect of the unconstrained pair-annihilation 
process even in the fast diffusion limit. We conclude that 
the system undergoes a mean-field transition, but with 
fluctuations in the particle number which deviate from 
those predicted by mean-field. 
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B. Spin-relaxation formulation 

We consider the system in the dual limit a — > oo 
where we can study the phase transition between the 
active phase and the absorbing phase in terms of the 
dimensionless variable u = (D — A)/(D + A). In the 
spin-relaxation picture this is the limit of fast Kawasaki 
spin-exchange where the system is spatially uncorrelated 
and hence completely characterized by the total magne- 
tization M = J2k °fc/2- The dynamics of the process (|l|) 
reduce in this limit to a random walk in the magnetiza- 
tion variable M, ranging from — L/2 to L/2. The master 
equation reads 

J t p M{t) = r M P M -i{t) + £m+iPm+i (t) 

-{r M +iM)PM{t). (31) 

The transition rates for this random walk with absorbing 
boundaries at M = ±L/2 are readily calculated as 

D + X ( 2Mu\ ( L 2 ,\ , x 

rM = 2L^{ 1 -—2){T- M ) 

D + X ( 2Mu\ ( L 2 ,\ , . 

iM = —2{ 1 + —2){-- M )- 

We find a bias towards to the boundaries M = ±L/2, 
i.e., the fully magnetized absorbing states with all spins 
up or down resp., for u < 0. For u > the system is 
biased to the center, corresponding to the active phase. 

For an initial state which is symmetric under spin-flip 
Si — > —Si the mean (M) vanishes for all times in both 
the active and absorbing phase and hence is not suit- 
able to characterize the system. For the same reason 
an naive mean-field approach by setting (erf (i)<x| (i) ) = 
( erf (t) ) ( er| (i) ) would not give any information on the 
dynamics of the spin-fluctuations. The quantity that 
characterizes the phase transition are the fluctuations 
(M 2 ) = ^2 M M 2 PM(t) in the magnetization, i.e. the 
mean-square displacement of the random walk. In the 
active regime this quantity is proportional to the system 
size whereas in an ordered state ( M 2 ) ~ L 2 . 

First consider the phase transition point u = 0. 
From the considerations above we know that the sta- 
tionary state is inactive. The only question of inter- 
est is the approach to stationarity from some random 
initial state. From (|l]) one obtains d/(dt)(M 2 ) = 
2A(L 2 /4 - ( M 2 ))/{L - 1) which is readily solved by 

(M 2 (t)) = L 2 /4 + ((Af 2 (0)) -X 2 /4)e- 2At/(L_1) . (34) 

The approach to the stationary value is exponential on a 
time scale 



For large system size and initial times t«L, the fluctu- 
ations in the magnetization grow linearly in time. 



For a / the equations of motion for the moments 
( M 2k ) are too complicated for direct analysis. We de- 
fine M = M/VL and study only the thermodynamic 
limit. Using the master equation ( |3l| ) the stationarity 
condition dj (dt)(M 2k ) = for the moments of M yields 
the recursion relation 

(M 2k ) = {2k- l)(Af 2fe - 2 )/(4u) (36) 

which shows that the stationary distribution in the active 
phase u > is Gaussian with variance l/(4u): 



P*(M) = 




This yields the final result 

(M 2 ) = l 1/(4U) U >n (38) 
N ' I oo u < v ' 

All other stationary moments in the active regime follow 
from the Gaussian nature (|37|) of the statistics. We read 
off a critical exponent k = 1 for the divergence of ( M 2 ) 
with u as the system approaches the critical point u = 0. 

IV. RELAXATIONAL BEHAVIOUR IN FINITE 
SYSTEMS 

The exact solution || for the dynamics of the spin-spin 
correlation function on the line A = D implies a crossover 
time t ~ L 2 from a power law relaxation to exponential 
relaxation to the absorbing state. One then expects this 
to hold throughout the inactive phase. 

On the other hand, any finite system has only one sta- 
tionary state, which is the empty lattice in particle lan- 
guage, corresponding to the magnetically ordered states 
with all spins up or all spins down. It is therefore of inter- 
est to study the relaxation towards this state in that re- 
gion of parameter space that constitutes the active phase 
of the infinite system. 

The precise location of the phase transition line is not 
known, but we know that the line A = (no unrestricted 
pair annihilation) belongs to the active phase and thus 
we may get some insight by studying the system in the 
immediate neighbourhood A <^ a, D of this line. To this 
end we adopt a similar strategy as in the previous sec- 
tion by assuming A to be so small that the system had 
sufficient time to relax to its A = stationary distri- 
bution between two successive pair-annihilation events. 
This limiting procedure can be made rigorous by taking 
a — ► oo and keeping A, D fixed. 

In this limit the system reduces effectively to a two- 
state system, i.e. the system is completely characterized 
by stating whether the system is empty (state | )) or not 
(which we denote by | 1 )). The latter state represents the 
stationary distribution of the system with A = in which, 
because of detailed balance for this reduced process, all 



G 



states with an even, non-zero number of particles have 
equal probability p — 1/[2^ L_1 ^ — 1]. The transition rates 
between these two states characterizing the system are 
then trivial to work out: The transition from | ) to | 1 ) 
is zero because | ) is an absorbing state. On the other 
hand, counting the number of states represented by 1 1 ) 
for which a pair annihilation event leads to the empty 
lattice yields a transition rate 



L/2-1 



l/r act = (AL)/^- 1 ) - 1] 



(39) 



for transitions from state 1 to state 0. Hence, at time 
t the system is in the absorbing state with probability 
Po(t) — 1 — e~ t / Tact and in each non-empty state with 
probability Pt(t) = e^l^ /p^" 1 ) - 1]. 

For the particle density this behaviour implies the ex- 
act result 



Pit) 



-t/T act 



(1/2Y 



(40) 



Because of the fast intermediate relaxation to the equi- 
librium state of the a = oo process the density and the 
density-correlations have no spatial dependence and the 
diffusion rate does not enter. The crossover time r ac t for 
reaching the absorbing state in the active region of the 
phase diagram for the infinite system (small A) diverges 
exponentially in system size. 

Finally, we study the dynamical behaviour of the sys- 
tem for large, but fi nite L in the limit of fast diffusion di- 
cussed in Sec. [II A . We recall relation (17) which relates 
the decay of the particle density to the survival probabil- 
ity of two neighouring particles in an empty lattice. This 
quantity can be interpreted as a first-passage-time distri- 
bution for two annihilating and branching random walk- 
ers: When two random walkers in an empty lattice an- 
nihilate for the first time, the dynamics stop. Therefore 
the density decay equals one half this first-passage-time 
distribution and the mean-first-passage-time (MFPT) 



lim (0\(H + c)- 1 \k,k+ 1) 

c — >0 



(41) 



gives the crossover time scale on which the system reaches 
the absorbing state. 

This quantity can be evaluated numerically for any 
point in parameter space by inverting the time evolution 
operator c + H for finite system size, then taking the ma- 
trix element (^l[) and finally calculating the limit c — > 0. 
For an analytical treatment for large D we note that the 
MFPT from some site k to an absorbing site k = for 
a general random walk with nearest neighbour hops on 
L + 1 sites can be expressed in terms of the hopping rates 
p3[ . In the mapping of the B ARW to the random walk 
the MFPT r is equal to the MFPT of the random walker 
starting at site 2. With the hopping rates (|26|), ( |27| ) we 
find after some rearrangement of terms 



XL ^ 



(42) 



fc=0 



with 



Cfe = 



L\k\T(\(L- 2)/a + l) 



(2k + 2)! (L — 2k — 2)\T(X(L - 2)/a + k + l) 



(43) 

We note first that in the limit a — > oo the MFPT co- 
incides with the relaxation time (|39|), as expected from 
duality. To study the asymptotic behaviour of r for finite 
a (active phase) we determine the value k for which cj, 
gives the largest contribution to sum on the r.h.s. of (|I|). 
We find k = p*L/2 with the stationary density given by 
(^8| ) . Using the Stirling formula for the Gamma- function 
and expanding Ck around fco yields for non- vanishing den- 
sity p* the asymptotic form of the crossover time 



(1 - 2 ( o*) (1-2p * )/(2p * ) 
(1 - p*)( 1 -p')/(p*) 



(44) 



up to subleading power-law corrections in system size. 
Therefore, in the active region of the phase diagram the 
crossover to absorption in a finite system takes place on 
a time scale which is exponentially large in system size, 
with a density-dependent amplitude. 

For p* = 0, i.e. in the absorbing phase , the MFPT can 
be read off directly from (|i^), (|43|), since only the term 
with k = contributes. One finds 



Tabs 



(L-l) 
2A 



(45) 



This power law differs from the crossover behaviour 
r abs ~ L 2 /D for finite diffusion constant D. The MFPT 
for this point in parameter space coincides with the re- 
laxation time (pq) in the dual point. 



V. FINAL REMARKS 



The duality relation (|13|) divides the parameter space 
into two distinct regions separated by the self-dual line 
( |l5| ) . Both regions are mapped onto each other and hence 
have the same physical properties. This is of practical 
usefulness for a numerical survey of the system and for 
the determination of the phase transition line since only 
part of the parameter space needs to be investigated. 

In particular, the line a = maps onto the line D = A 
and the fast-diffusion limit to the limit a — ► oo. The ob- 
servation that for large D any small a brings the system 
into the active phase translates into a phase transition at 
D = A in the limit a — > oo. This is a new result which 
definitely clarifies the unresolved issue of the location of 
the phase transition line for large a. Numerical analy- 
sis of the model for large rates a is reported to be very 
difficult PJ2C|]. Our exact result confirms the conjecture 
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of Ref. H on the location of the phase transition point 
in this numerically untractable limit. In the 5 — p ex 
phase diagram of Ref. j^] the limit a — > oo corresponds 
to p ex — > 1 and the phase transition point D — X cor- 
responds to 6 = 0. The dual limit D — > oo covers the 
neighbourhood of the point S = —l,p ex = 0. Our result 
translates into an infinite slope of the phase transition 
line in this representation at this point. In the active 
phase of this region the exact stationary particle num- 
ber distribution is Gaussian with a stationary density p* 
given by the mean- field value (p3f). The density fluctu- 
ations ([BO]) deviate from the mean-field result (|24|) by a 
factor 2(1 - p*). 

On the self-dual line we find from ( |l7| ) that the density 
expectation value pk it) equals one-half the survival prob- 
ability at time t of two particles placed initially at two 
neighbouring sites in an otherwise empty lattice. Hence 
the phase transition from the absorbing phase to the ac- 
tive phase may be rephrased as a mean-first-passage-time 
(MFPT) transition for random walkers which branch and 
annihilate. We expect the MFPT in a finite system to 
change from a power-law divergence (in system size) to 
an exponential divergence not only at infinite diffusion 
rate (Sec. IV), but also for finite D. Thus this numer- 



ically accessible quantity provides an alternative way of 
determining and characterizing the PC phase transition. 



VI. ACKNOWLEDGEMENTS 

We thank J. Cardy, N. Menyhard, Z. Racz and U. 
Tauber for useful discussions. G.M.S. thanks the Insti- 
tute for Theoretical Physics of the Eotvos University for 
hospitality and partial support (Grant OTKA OTKA T 
019451). 



[8] M. Droz, Z. Racz and J. Schmidt, Phys. Rev. A 39, 2141 
(1989). 

[9] B. P. Lee, J. Phys. A 27, 2633 (1994). 
[10] E. Siggia, Phys. Rev. B 16, 2319 (1977); S. Sandow and 
S. Trimper, Europhys. Lett. 21, 799 (1993); F. C. Al- 
caraz, M. Droz, M. Henkel and V. Rittenberg, Ann. Phys. 
(N.Y.) 230, 250 (1994); G. M. Schiitz, J. Stat. Phys. 79, 
243 (1995). 

[11] G. M. Schiitz, Integrable Reaction- Diffusion Processes 
and Quantum Spin Chains (to be published). 

[12] F. Spitzer, Adv. Math. 5, 246 (1970). 

[13] J. W. Evans, in Nonequilibnum Statistical Mechanics in 
One Dimension, ed. V. Privman (Cambrige University 
Press, Cambridge UK, 1996). 

[14] D. C. Torney and H. M. McConnell, J. Phys. Chem. 87, 
1941 (1983). 

[15] J. E. Santos J. Phys. A 30, 3249 (1997). 

[16] D. Levy, Phys. Rev. Lett. 67, 1971 (1991). 

[17] G. Schiitz J. Phys. A 26, 4555 (1993). 

[18] G. M. Schiitz, Z. Phys. B 104, 583 (1997); G. M. Schiitz 
and K. Mussawisade, Phys. Rev. E 57 (1998) (in press). 

[19] Such a relationship between a stochastic processes is 
called an enantiodromy relation, as opposed to a similar- 
ity transformation like the domain- wall duality which re- 
lates one stochastic Hamiltonian to another, rather than 
to its transposed. 

[20] N. Menyhard and G. Odor J. Phys. A 28, 4505 (1995). 

[21] D. Kandel, E. Domany and B. Nienhuis, J. Phys. A 23, 
L755 (1990). 

[22] This and the other limits of fast rates discussed in this pa- 
per can be treated rigorously by taking the limit of large 
rates in the formal solution of the master equation, see 
for details. 

[23] G. H. Weiss, J. Stat. Phys. 24, 587 (1981); K. P. N. 
Murthy and K. W. Kehr, Phys. Rev. A 40, 2082 (1989), 
and erratum in Phys. Rev. A 41, 1160 (1990). 



[1] T. E. Harris, Ann. Prob. 2, 969 (1974); H. Takayasu and 
A. Yu. Tretyakov 68, 3060 (1992); I. Jensen, Phys. Rev 
E 47, Rl (1993), Phys. Rev. Lett. 70, 1465 (1993). 

[2] P. Grassberger, F. Krause and T. von der Twer, J. Phys. 
A 17, L105 (1984); P. Grassberger, J. Phys. A 22, L1103 
(1989); M. H. Kim and H. Park, Phys. Rev. Lett. 73, 
2579 (1994); I. Jensen, J. Phys. A 26, 3921 (1993), Phys. 
Rev. E 50, 3623 (1994); D. ben-Avraham, F. Leyvraz 
and S. Redner, Phys. Rev. E 50, 1843 (1994). Some of 
the references consider BARW's with both even and odd 
number of offsprings. 

[3] N. Menyhard, J. Phys. A 27, 6139 (1994). 

[4] J. Cardy and U. Tauber, Phys. Rev. Lett. 77, 4780 
(1996); J. Stat. Phys. 90 (1998) (in press). 

[5] R. J. Glauber, Math. Phys. 4, 294 (1963). 

[6] K. Kawasaki, Phys. Rev. 145, 224 (1966). 

[7] Z. Racz, Phys. Rev. Lett. 55, 1707 (1985). 



8 



